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ABSTRACT 

We present a catalog of 49 spectroscopic strong gravitational lens candidates selected from a Sloan 
Digital Sky Survey sample of 50996 luminous red galaxies. Potentially lensed star-forming galaxies 
are detected through the presence of background oxygen and hydrogen nebular emission lines in 
the spectra of these massive foreground galaxies. This multiline selection eliminates the ambiguity of 
single-line identification and provides a very promising sample of candidate galaxy-galaxy lens systems 
at low to intermediate redshift, with foreground redshifts ranging from 0.16 to 0.49 and background 
redshifts from 0.25 to 0.81. Any lenses confirmed within our sample would be important new probes 
of early-type galaxy mass distributions, providing complementary constraints to those obtained from 
currently known lensed high-redshift quasars. 

Subject headings: gravitational lensing — galaxies: elliptical and lenticular, cD — galaxies: starburst 
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1. INTRODUCTION AND MOTIVATION 

In its strongest form, gravitational lensing produces 
unmistakably distorted, amplified, and multiplied images 
of distant astronomical objects. It is therefore not sur- 
prising that the majority of known galaxy-scale gravita- 
tional lens systems have been discovered through imag- 
ing observations. However, a small number of lenses 
have been discovered spectroscopically, with the spec- 
trum of a targeted galaxy showing evidence of emission 
from a background source and follow-up imaging reveal- 
ing lensing morphology. The three most secure exam- 
ples are the lensed quasars 223 7+0305 llHuc hra et ah! 
M98l and SDSS J0903+5028 (|.Tormston et a,U 12003) 
and the lensed Lvman-a emitt ing galaxy 0047—2808 
l|Warren et all 119961 Il998l 119991). Several authors have 
made predictions for the frequency of lensed quasar 
discoveries in galaxy redshift surveys ijKochanekl 119921: 
IMortlock fc Wer^l2T)00l Eool. Others, inspired by 
the discovery of 0047—2808, have undertaken spec- 
troscopic searches for lensed Lyman-q-bright galaxies 
llHewett et all 1200(1 iWillis et alJl2007i lHall et id] 1200(1 
iBurles et alJ l2000f) . The idea behind such searches is 
that a massive foreground galaxy should act as an ef- 
fective gravitational lens of any object positioned suffi- 
ciently far behind it and at small enough impact param- 
eter, and any emission features from such lensed objects 
should be detectable in the spectra of the foreground 
galaxy. Therefore a search for discrepant emission fea- 
tures in galaxy spectra can potentially lead to a sample 
of gravitational lens systems that would not be discov- 
ered in broadband imaging searches due to faintness of 
source relative to lens. This "lenses-looking-for-sources" 
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approach is complementary to lens searches such as the 
recently completed Cosmic Lens All-S ky Survey (CLASS; 
iMvers et alJl2003l iBrown et alJ t2003) which proceed by 
targeting sources and looking for evidence of an inter- 
vening lens. 

With its massive scale and quality of data, the spec- 
troscopic component of the Sloan Digital Sky Sur- 
vey (SDSS) provides an unprecedented opportunity for 
spectroscopic galaxy-galaxy gravitational lens discovery. 
This paper presents the first results of such a search 
within a sample of ~ 5 1,0 00 SDS S Lum inous Red Galaxy 
(LRG) spectra l|Eisenstein et aLlfeOOll hereafter E01): a 
catalog of candidate lensed star-forming galaxies at in- 
termediate redshift. The lens candidates presented in 
this paper were detected by the presence of not one but 
(at least) three emission lines in the LRG spectra iden- 
tified as nebular emission from a single background red- 
shift: [O n] AA3728 and two out of the three of H/3 A4863, 
[O III] A4960, and [O m] A5008 (all emission-line wave- 
lengths in this paper are vacuum values, in Angstroms). 
This implies a maximum redshift of z ~ 0.8 for any candi- 
date lensed galaxies: at higher redshifts, [O in] emission 
moves off the red end of the SDSS spectrograph. 

The lensing cross section of a particular foreground 
galaxy is lower for intermediate-redshift sources than 
high-redshift sources, and in this sense the lens search 
we describe here is at a disadvantage relative to searches 
for lensed high-redshift Lyman-a emitters. However, 
the identities of the emission lines in the sample we 
present here are absolutely unambiguous. This can- 
not be said for spectroscopic Lyman-a lens candidates, 
which are typically detected as single discrepant emission 
lines and are difficult to distinguish from lower-redshift 
emission or (in the case of a huge survey such as the 
SDSS) non-astrophysical spectral artifacts. (A forth- 
coming paper will present SDSS single-line [O il] and 
Lyman-a lens candidates.) Our spectroscopic lens search 
based on multiple-line detection has further advantages. 
First, the source redshift of any lensed galaxies will be 
known from the outset. Second, for a given limiting 
line flux, star-forming galaxies at intermediate redshift 
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are more numerous on the sky than Lyman-g emitters 
at hi gh redshift (see iHippelein et alJ l200.1t iMaier et al.1 
2003) , and hence strong lensing events could be more fre- 
quent despite decreased lensing cross sections. The non- 
lensed source population shoul d thus be more am enable 
to study and characterization (Hogg et alJll998l for ex- 
ample), facilitating lens statistical analysis. Finally, 
intermediate-redshift lensed galaxies would probe the 
mass distribution of the lens population in a systemati- 
cally different manner than do high-redshift sources. 

The search for new lens systems is motivated by 
the unique power of gravitational lensing to con- 
strain the mass distrib u tion o f the lensing galaxy. 
iMiralda-Escude fe Leharl l)1992|) argued that strong 
galaxy-galaxy gravitational lenses in the optical band 
should be both ubiquitous and of great scientific inter- 
est, but their detection has proved difficult. In addi- 
tion to the optical E i nstein ring 0047—2808 mentioned 
above, IHewett et al.l l|2000f) have published one more 
spectroscopic galaxy-lens candidate. Hubble Space Tele- 
scope (HST) imaging has also had so me success in de- 
tectin g lensed galaxies, as described bvlCrampton et al.1 
112001. iR^tunga. Griffiths, fc OstrandeH l|1999D. and 
iFassnacht et all ()2003|) . but spectroscopic confirmation 
remains a challenge in most cases. Lensed galaxies are 
especially useful for investigating the mass distribution of 
the lensing galaxy because unlike lensed quasars they are 
spatially res olved sources which provide more detailed 
constraints l|Kochanek fc Naravanl IT992. for example). 
As with lensed quasars, the relative angular deflections 
between multiple images can constrain the total deflect- 
ing lens mass. But the relative distortions of multiple 
images can provide constraints as well, probing the slope 
of the radial mass profile and allowing more detailed con- 
nection with theories of galaxy formation and Cold Dark 
Matter (CDM). The promise of our candidate sample is 
enhanced by the fact that lens and source redshifts for 
any confirmed lens systems will be known from the out- 
set with great precision; this is generally not the case for 
lenses discovered in optical imaging or radio data. 

The lack of sufficient lensing constraints on the radial 
slope of the galaxy mass profile in lensed quasar sys- 
tems is the greatest source of uncertainty in the Hubble 
constant (Hp) values derived from the time delays be- 
tween images l|Schechterll2004T) . which remain systemat- 
ically low compared to the current best Hp values from 
the distance ladder llFreedman et alJ| 2001l) a nd cosmic 
microwave back ground IjSoergel et al.ll2003[k iKochanekl 
(2002a b, 2003) highlights the importance of resolving 
this discrepancy. Any galaxy-galaxy lenses confirmed 
within our current list of candidates would help to solve 
this problem by expanding the sample for ensemble lens- 
ing studies of early-type galaxy mas s profiles like that of 
IRusin. Kochanek. fc Keetonl (|2003bT) and allowing for the 
inclusion of distortion information as mentioned above. 
The generally lower lens and source redshifts of our can- 
didate systems as compared to known lens systems are 
also of particular interest because the images of any 
strong lenses within our sample will tend to form at 
smaller physical radii wi thin the lens galax y than in 
current lens systems (see IRusin et al J r2003al Table I), 
thus probing elliptical galaxy mass profiles in a system- 



atically different manner 6 . A large enough sample of 
strongly lensed star-forming galaxies could even be mon- 
itored directly for supernovae. Multiple images of the 
same supernova would give time delays and Hp mea- 
sure ments, as numerous a uthors have envis ioned (pj efsdal 
T9fi ^|Porcia,ni fc Madaull27)00t IhT^ EtoI IGooharTIHTI 
2002: IQguri fc Kawand 12003^ 7 If we assume one super- 
nova per galaxy per century as a n order-of-magnitude 
estimate (Doggett fc Brancblll985fl . a program to moni- 
tor ~ 50 lensed galaxies for several years would have a 
good chance of success. 

Numerous studies have addressed the inability of 
standard macroscopic galaxy mass models to pre- 
dict the flux ratios observed between images of 
strongly lensed quasars. Some authors have suggested 
that this phenomenon may be due to the effects 
of substructure in the lensi ng galaxy on the scale 
of either CDM sub-halos jMao fc Schneiderl Il998t 
Metcalf fc MadavJl200l| lChiball2002t iDalal fc Kochanekl 
2002ft or stars Jwitt. Mao, fc Schechterl 119951 
Schechter fc Wa mbsganss 2002), while others have 
invoked macroscopically mo dified lens models 
as a possible expla n ation l|Evans fc Witti l2003t 
iMoller. Hewett. fc Blainl 12003ft . A sizeable sample 
of strongly lensed galaxies would provide valuable 
insight into this problem: by virtue of the arcsecond- 
scale intrinsic angular size of the source, galaxy-galaxy 
lenses will be immune to the lensing effects of any sub- 
structure with lensing deflection angles on the scale of 
microarcseconds (stars) or milliarcseconds (sub-halos). 
If the suppressed saddle point image fluxes seen in 
many quasar lenses are indeed due to macroscopic 
inadequacies in the standard lens models, similar effects 
should be seen in galaxy-galaxy lenses. If on the other 
hand these effects are due to micro/millilensing by 
substructure, they should be absent from lensed galaxy 
systems. 

2. SEARCH SAMPLE 

The Sloan Digital Sky Survey is a project to image 
roughly one-quarter of the sky in five optical bands 
and obtain spectroscopic follo w-up observ a tions of ~ 
10 6 galaxies and 10 5 quasars. lYork et aD JgQflCj ft pro- 
vide a technical summary of t he survey. iGu nn el; all 
(11998ft desc ri be the SDSS c amera. iFukugita et al.l iTl996). 
iHogg et all 1)200 1|) . and ISmith et alJ (12002ft discuss 
the p hotometric system and c alibration, [ pier eij alJ 
( 2003) discuss SDSS astrometry, iBlanton et all l|2003ft 
present the spe ctroscopic plate tiling algorithm, and 
iStoughton et alJ l|2002ft and lAbazaiian et alJ l|2003ft de- 
scribe the survey data products. Approximately 12% 
of the galaxy spectroscopic fibers are allocated to the 
LRG sample (E01), selected to consist of very lumi- 
nous (> 3L*) — and hence massive — early-type galaxies 
at higher redshift th an the galaxies of the main sample 
(jStrauss et al.ll200'2^ . We expect these LRGs to be par- 
ticularly effective gravitational lenses of any objects po- 
sitioned suitably behind them, and we concentrate our 
spectroscopic lens search on them. We also note that 
LRGs should have little dust and therefore any lensed 
background galaxies should suffer minimal extinction. 

6 The lowest-redshift l ensed quasar cu rrently known, RXJ 
1131-1231, has z = 0.658 ISluse et alj|2003T) . 
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The sample for our current study consists of 50996 spec- 
tra taken between 5 March 2000 and 27 May 2003 of 
SDSS imaging obj ects flagged as GALA XY_RED by the pho- 
tometric pipeline l)Lupton et al.ll2001[) for passing either 
of the LRG cuts described by E01, reduced by the SDSS 
spectroscopic pipeline (J. Frieman et al., in preparation), 
and selected to have redshifts between 0.15 and 0.65 
as determined by the specBS rcdshift-finding software 
(D. J. Schlegel et al., in preparation). The low-rcdshift 
cutoff is needed because less massive galaxies start to 
pass the photometric cuts below z < 0.15 and pollute 
the volume-limited LRG sample; see the discussion in 
E01 and the bimodal L RG sample redshift histogram in 
iStoughton et all (|2002L Fig. 14). 

In addition to being much more massive than the av- 
erage galaxy, LRGs have another property that makes 
them well suited to a spectroscopic lens survey: their 
spectra are extremely regular and well-characterized (see 
lEisenstein et al J feOOffl . To determine the spectroscopic 
redshift of an SDSS target galaxy with observed specific 
flux f\ and one-sigma sky+source noise spectrum <j\ , the 
specBS program employs a small set of galaxy eigenspec- 
tra (four in the reductions for this study) derived from 
a rest-frame principal-component analysis (PCA) of 480 
galaxy spectra taken on SDSS plate 306, MJD 51690. 
This eigenbasis is incrementally redshifted, and a model 
spectrum is generated from the best-fit linear combina- 
tion to the observed spectrum at each trial redshift, with 
the final redshift assignment given by the trial value that 
yields the overall minimum % 2 . Although redshift is the 
primary output of this procedure, a byproduct of specBS 

is the best-fit model spectrum itself, fx- In the case of 
LRGs, f\ typically provides a very detailed and accurate 
fit to f\, with a reduced \ 2 of order unity over almost 
4000 spectral pixels (roughly 1600 spectral resolution el- 
ements) attained with only 8 free parameters: a redshift, 
the four eigen-galaxy coefficients, and the three terms 
of a quadratic polynomial to fit out spectrophotomet- 
ric errors and extinction effects, both of which exist at 
the few-percent level 7 . This extremely regular spectral 
behavior allows us to form residual LRG spectra 

= fx - fx (1) 

that are in principle realizations of ax- Nebular emis- 
sion lines from galaxies along the line of sight other than 
the target LRG will not be modeled by specBS and 
should appear as significant features localized in wave- 
length within these residual spectra. Figure Q] shows the 
median l-a line flux sensitivity within our LRG residual 
spectrum sample as a function of wavelength. The (20th, 
50th, 80th)-percentile LRG spectra themselves have a 
median signal-to-noise per pixel of (3.3, 5.1, 9.6) at the 
SDSS resolution of A/AA « 1800. 

3. CANDIDATE SELECTION 

This section describes our candidate selection routine 
in detail. Briefly stated, we select as initial candidates 
those spectra that show both blended [O n] AA3728 at 
S/N > 3 and two out of the three lines H/3 A4863, [O III] 
A4960, and [O in] A5008 at S/N > 2.5, then cull the 

7 The specBS outputs for the SDSS-DRl sample are available 
from the website http://spectro.princeton.edu 
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Fig. 1. — LRG sample median line-flux noise spectrum. Shown 
is the l-a noise on best-fit line fluxes for optimal extraction of 
Gaussian-shaped residual emission features with spectral width 
cr =1. 2-pixel (~ 83 km s" 1 ). Reported pixel flux variances have 
been rescaled as described in Appendix [A] prior to the calculation 
of this noise spectrum. 



candidate list by applying cuts based on more detailed 
fits to the presumed emission features, and finally remove 
any obviously spurious detections. This selection process 
yields a substantial number of promising systems without 
an excess of obvious false positives. 

3.1. Initial Emission Feature Detection 

The key element in the first step of our lens candidate 
selection (described fully in § 13. 21 below ) is a straightfor- 
ward matched-filtering procedure to search for significant 
emission features by fitting a G aussian line profile at each 
point in the residual spectrum l)Prattll978HHewett et alJ 
I1985|) . We describe our implementation here so as to be 

explicit. Let /j 1 *"* be the residual flux in pixel j and <r| be 

the statistical variance of /j r) 8 . Also let {ui\ describe 
a Gaussian kernel, centered on i = 0, with i running 
from —i]i m to i]i m , and normalized such that Ui = 1. 
The maximum-likelihood estimator Aj for the line flux 
Aj of any {wi}-shaped residual emission feature centered 
on pixel j is that which minimizes 

x?= £ (^-/Sii))7-u • (2) 

i=— ' 

Differentiating J2J with respect to Aj , setting the result- 
ing expression to zero, and solving yields 

A,=cf , (3) 

where we have defined the convolutions 

Cf >s E /SV/^-H) ■ ( 4 ) 

i 
i 

8 Conversion from units of ergs cm -2 s -1 A -1 to units of ergs 
cm -2 s — 1 pixel -1 is made using the re-binned SDSS spectroscopic 
pixel scale relation d\ = A X 10~ 4 ln(10) d(pixcls). 
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The variance of Aj (under the assumption of uncorre- 
cted Gaussian errors in as described by crj) is given 

by 

°h = ^/ c ? ] • (6) 

The signal-to-noise ratio for a fitted Gaussian profile cen- 
tered on pixel j is therefore 

(5/JV)j = Cf ) /yftif* ■ (7) 

Our null hypothesis is an absence of emission features 
in the residual spectra that should manifest as the 
{{S/N)j} being Gaussian-distributed about zero with 
unit variance: this should hold at most wavelengths in 
most spectra. We approach the initial search for emission 
lines in the residual spectra as a search for significance 
peaks with (S/N) greater than some threshold value. Al- 
though insensitive to goodness-of-fit, this convolution- 
based detection scheme executes quickly and is therefore 
well suited to the initial search for residual emission fea- 
tures within our large spectral sample. Appendix ^de- 
scribes a noise-rescaling process that we employ to con- 
trol the incidence of false-positive emission- feature detec- 
tions (due primarily to imperfect sky subtraction) with- 
out masking regions of the spectrum. 

3.2. Multi-Line Background Systems: Detection, 
Fitting, and Rejection 

Multiple emission features at the same redshift will 
have redshift-independent wavelength ratios. The fully 
reduced SDSS spectra have been re-binned at a constant- 
velocity pixel scale of 69 km s _1 , giving a redshift- 
independent pixel offset between features. Our oper- 
ational scheme is thus to search for coincident (S/N) 
peaks between multiple copies of a single filtered resid- 
ual spectrum that have been shifted relative to one an- 
other. For [O n] AA3728 detection, we filter each residual 
spectrum with a a = 2.4-pixel Gaussian kernel (matched 
to the typical width of blended [O n] AA3728 emission 
seen in SDSS starburst galaxies). We take copies of the 
same residual spectrum "blueshifted" by integer pixel 
amounts so as to place RS A4103, H 7 A4342, H/3 A4863, 
[O in] A4960, [O in] A5008, [N n] A6550, Ha A6565, [N n] 
A6585, [S n] A6718, and [S n] A6733 as close as possible to 
3728.48 A, the geometric-mean wavelength of the [O II] 
AA3728 doublet. These shifted spectra are filtered with 
a a = 1.2-pixel Gaussian kernel (matched to the typical 
width of SDSS starburst [O III] A5008 emission), with 
the sub-pixel part of the line offset relative to 3728.48 
A incorporated by offsetting the kernel. Any pixel in 
the filtered S/N spectra with value greater than 3 for 
[O n] AA3728 and value greater than 2.5 for two out of 
H/3 A4863, [O III] A4960, and [O Hi] A5008 is tagged as 
a "hit". A group of adjacent "hit" pixels is reduced to 
the single pixel with the greatest quadrature-sum S/N 
for lines detected above the threshold (in effect, the pixel 
most inconsistent with the null hypothesis). Spectra with 
more than one isolated hit are rejected. The spectra are 
only searched in regions that would correspond to emis- 
sion from > 5000 km s _1 behind the targeted LRG. 

The preceding selection leads to 163 hits within our 
51,000 spectra. For each hit, we explore a grid of redshift 



and intrinsic emission-line-width values for the back- 
ground galaxy to find a best-fit model. At each grid 
point we fit a Gaussian profile to any emission line ini- 
tially detected above a 2.0-S/N threshold, with the line 
center determined by the trial redshift and line-width 
given by the quadrature sum of the trial intrinsic line- 
width and the measured spectrograph resolution at the 
observed wavelength of each line. [O n] AA3728 is fit 
with a double Gaussian profile. We adopt as best values 
for background redshift zbg and (Gaussian-c) intrinsic 
line-width o~ e those that give the minimum x 2 over all 
detected lines. The Zbg extent of our grid corresponds 
to ±2 pixels, and the explored o~ e range runs from to 2 
pixels (0 to 138 km s -1 ). 

Following these fits, we subject the candidate sample 
to several cuts. First, we reject any system where no 
convergent (minimum x 2 ) value for zbg is found within 
the explored ±2-pixel range: real background emission 
systems should have convergent best-fit redshifts roughly 
coincident with their convolution-based S/N peaks. Sim- 
ilarly, we cut systems with no convergent o~ e between 
and 2 pixels: the features of real systems should neither 
be narrower than the spectrograph resolution nor wider 
than a reasonable upper limit for emission features of 
faint star-forming galaxies. Finally, we compute a to- 
tal signal-to-noise ratio for the fit, defined as the total 
best-fit flux in [O n] AA3728 and all other lines initially 
detected at S/N > 2.5 divided by the quadrature-sum 
of the T<7 noise from those line fits, and impose a cut 
in the total-5/iV-Xr plane (x 2 being the x 2 P er degree 
of freedom in the fit). We cut any system with a total 
S/N less than the greater of 6 and 6 + 3(%^ — 1). This 
removes both low-S/N candidates and candidates whose 
X 2 values are too high to be explained by high-S/N emis- 
sion features showing significant non-Gaussian structure. 
(These cuts were designed to be a quantitative expression 
of our own judgements about which candidate systems 
look real upon spectrum inspection and which do not; 
the justifications are somewhat a posteriori.) This cut- 
ting procedure reduces the 163 hits to 61 candidate sys- 
tems. Finally, we prune 12 candidates from the list that 
survive the automated culling but are clearly explained 
by either over-fit LRG stellar absorption, under-modeled 
LRG line emission, exceptionally poor data quality, or a 
generally flawed template fit, leaving the 49 good candi- 
date systems that we present. 

This search for background galaxy emission lines digs 
rather deep into the noise of our spectroscopic sample. 
To gauge the incidence of false positives in our final can- 
didate list, we make a parallel run of the detection, fit- 
ting, and automated rejection procedure with the follow- 
ing rest-wavelength perturbations: H(3 A4863 — > 4833, 
[O m] A4960 -> 4945, and [O ill] A5008 -> 5023. These 
perturbations alter all of the redshift-independent wave- 
length ratios among these lines and between all of them 
and [O n] AA3728; this modified detection procedure no 
longer selects for real multi-line emission, but only for 
noise features. The "false candidates" that result from 
this perturbed procedure are randomly shuffled along 
with the candidates from the original procedure, and all 
are examined together when making the final pruning 
judgements. The perturbed procedure yields 88 hits and 
7 post-cut candidates; all 7 are pruned upon inspection 
without knowledge of their intrinsic falseness. This im- 
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plies that the vast majority of our candidates are indeed 
background galaxies and not simply noise features. 

4. CANDIDATE SYSTEMS 
4.1. Catalog 

In this section we present our catalog of 49 candi- 
date lensed star-forming galaxies selected to have [O II] 
AA3728 emission at S/N of 3 or higher and emission 
from two out of the three of H/3 A4863, [O ill] A4960, 
and [O in] A5008 each at S/N of 2.5 or higher, at a 
rcdshift significantly greater than that of the primary 
target LRG. Table lists various properties of the can- 
didate lens systems. LRG r de Vaucouleurs model mag- 
nitudes are determined from SDSS imaging and photo- 
metric reduction. LRG redshifts and velocity disper- 
sions are as provided by specBS; the software fits for 
velocity dispersions a v using a set of 24 stellar eigen- 
spectra derived from a PCA of th e ELODIE spectral 
library ijPrueniel fc Soubiranll200ij) . We report all a v 
values from the database, although some are likely unre- 
liable; see the notes of Tabled We also report emission- 
line redshifts of the detected background galaxies. Us- 
ing the observed LRG and background redshifts and 
the observed LRG <r v , and assuming a singular isother- 
mal sphere (SIS) LRG luminous+dark matter distribu- 
tion, we calculate a "best guess" for the angular scale 
of any lensing that might be present in these systems as 
A9 = 8ir(aZ/c 2 ){D L s/Ds). {D LS and D s are angular- 
diameter distances from lens to source and from observer 
to source 9 ). This is the separation between the two im- 
ages of a strongly lensed object in the SIS model; it is 
also the radius of the strong-lensing region of the image 
plane, and twice the radius of ring images of compact 
sourc es directly behind the lens l|Naravan &: Bartelmannl 
119961 for example). Finally, for each candidate system 
we report the detected background [O n] AA3728 line flux 
from the best-fit double Gaussian profile, a list of other 
lines detected at S/N > 2.5, and the total S/N of the 
background galaxy emission-line fit as defined in § 13.21 
The reported background line fluxes should be viewed 
with suspicion: the spatial alignments of the background 
galaxies are unknown, and the spectroscopic fibers will 
not necessarily have captured all of their light. In addi- 
tion, most of the reduced candidate LRG spectra have 
been rescaled based on smear exposures that attempt 
to account for LRG flux falling outside the 3" fiber 
l)Stoughton et al.ll2 002) ; we have removed this correction 
from the tabulated line fluxes, as it has no physical rela- 
tion to the background emission-line galaxies. Figure |21 
shows the SDSS discovery spectra and best-fit model 
spectra, along with close-up views of the residual (data 
— model) spectra in the wavelength ranges correspond- 
ing to redshifted background [O n] AA3728, H/3 A4863, 
[O in] A4960, and [O in] A5008. Smear corrections have 
not been removed from these plots. 

Although we detect line emission clearly, evidence of 
background galaxy continuum in the residual spectra of 
our candidate systems is scarce. This is not surprising, 
for three reasons. One, the LRG sample was selected 
for particular broadband color and luminosity, and sig- 

9 Throughout this paper, we assume a cosmology of (Qm , Oa) = 
(0.3,0.7) and a Hubble constant Ho of 70 km s — 1 Mpc -1 . Since 
A9 is proportional to a distance ratio, it is independent of Hq. 



nificant background continuum would likely perturb an 
LRG out of the sample. Two, any faint background con- 
tinuum present in an LRG spectrum will largely project 
onto the LRG-redshift eigenspectrum set and low-order 
polynomial fit used by specBS, and will be subtracted 
along with the LRG model when forming the residual 
spectrum. Three, these background galaxies are likely to 
be high-equivalent-width star- forming systems, and since 
their line fluxes are detected just above the noise thresh- 
old, the associated continuum will typically be lost in 
the noise. Nevertheless, we may obtain a higher signal- 
to-noise picture of the background galaxies that we de- 
tect by constructing a median spectrum as follows. First 
we transform the residual spectra of our candidate sys- 
tems into units of erg cm~ 2 s _1 pixel -1 , which is a 
rcdshift-indcpcndcnt quantity since the SDSS pixels are 
of constant velocity width. We then shift these resid- 
ual spectra into the rest frame of the background galaxy, 
rounded to the nearest whole pixel, and transform back 
to erg cm' - 2 s -1 A -1 . Next we renormalize the spectra by 
dividing each one by its best-fit [O ii] AA3728-flux value. 
We then take the median value at each pixel, and restore 
physical normalization by multiplying this median spec- 
trum by the sample-median best-fit [O n] AA3728-flux 
value. The resulting median spectrum is shown in Fig- 
ure 01 Although there is no discernible continuum in the 
individual residual spectra, we can see a 4000-A contin- 
uum break in the median spectrum; we also see H7, Ha, 
[N 11] and [S 11] emission lines in addition to the lines for 
which we select. This gives further evidence that we have 
successfully detected and identified background emission 
features. 

4.2. Lenses or Not? 

We have with certainty detected emission from galaxies 
behind foreground LRGs. For a system to be a strong 
gravitational lens, the background galaxy must be lo- 
cated at sufficiently small impact parameter relative to 
the LRG center. The true incidence of lensing within our 
sample can best be determined and studied with either 
narrow-band imaging or integral-field spectroscopy. Such 
observations could spectrally isolate the background line 
flux and resolve it spatially to reveal any lensing mor- 
phology. SDSS broadband imaging offers some hope for 
answering the lensing question, but in general Sloan im- 
ages do not detect late- type galaxies at redshift z ~ 0.5 
at very high S/N. Figure H shows 40" x 40" SDSS r 
postage-stamp images centered on the candidate sys- 
tems, with linear gray-scaling from — 3-er sky noise to 
peak LRG surface brightness. Evidence of significantly 
offset neighboring broadband emission is seen in some 
images, but it would be difficult to rule out many sys- 
tems as definite non-lenses based on SDSS-quality im- 
ages. In the spirit of a purely spectroscopic survey, we 
present as candidates all systems selected spectroscopi- 
cally. Furthermore, we note that many of these LRGs live 
in high-density group/cluster environments, and neigh- 
boring images may be at the LRG redshift and not the 
source of the background emission that we detect. 

In preparation for integral-field spectroscopic follow- 
up, we have obtained broadband reconnaissance images 
of 14 of our candidate systems (as indicated in Tabled 
using the Magellan consortium's 6.5-m Clay telescope 
at Las Campanas Observatory. These short g- and i- 
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Table 1. Properties of Candidate Lens Systems 
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a Names give truncated J2000 RA and Dec in the format HHMMSS.ss+DDMMSS.s. 

b A8 column gives an estimated angular lensing scale for each system, calculated as described in the text. 
c [O n] AA3728 line flux in units of 10 -17 ergs cm -2 s _1 . 

d Other lines detected at 2.5-cr or greater. Letters specify lines as follows: a = H5 A4103, b = H7 A4342, c = H/3 A4863, d = [O III] A4960, 
e = [O hi] A5008, f = [N 11] A6550, g = Ha A6565, h = [N 11] A6585, i = [S 11] A6718, j = [S 11] A6733. 
e See !j |4.2l for remarks on follow-up imaging of these systems. 

f Velocity dispersion (cr v ) value possibly unreliable: either too high (> 420 km s — 1 ), too low to be properly resolved (< 100 km s — 1 ), or 
LRG spectrum median signal-to- noise per pixel too low (< 10). 

g SDSS-DR1 contains lower signal-to-noise spectra of these systems from earlier observations, with Plate-MJD-Fiber values of 414-51901- 
266 (for SDSS J0324) and 612-52079-607 (for SDSS J1521). 
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Fig. 3. — Median rest-frame spectrum of background galaxies 
detected in our residual LRG spectra, created as described in the 
text and smoothed with a 5-pixel boxcar. Note the continuum 
break at 4000A and the absorption feature at 5270A, which are 
not discernible in the individual residual spectra. Also evident are 
the emission lines H7 A4342, Ha A6565, [N n] A6550, [N 11] A6585, 
[S 11] A6718, and [S 11] A6733, in addition to the emission lines for 
which we select. 



band exposures (120s each) were taken on the nights of 
2003 August 1 and 2 with the Magellan Instant Cam- 
era (MagIC) CCD im ager and active telescope optics 
ijSchechter et alJl2003|) . Conditions were similar on both 
nights: non-photometric due to cirrus clouds, and with 
~ 0'.'8 median FWHM seeing. For the majority of these 
systems, the images present no significant evidence for 
or against the presence of strong lensing; we describe the 
exceptions here. One system, SDSSJ0035, is almost cer- 
tainly not a strong lens: a relatively <?-bright galaxy can 
be seen clearly in the images, approximately 2'.'6 offset 
from the LRG. This far exceeds the anticipated lensing 
scale of about 0'.'4 reported in Table Q The images of 
SDSSJ2147 reveal a bluer galaxy approximately 1'.'3 to 
the South, about at the limit of the estimated lensing 
scale. We see no particular evidence for strong lensing, 
but deeper observation of this system could be of in- 
terest. Two systems, SDSSJ0037 and SDSSJ0216, show 
evidence of multiple imaging. Figure [S] shows <?-band 
images of these two systems, along with difference im- 
ages created by smoothing aligned z-band images with a 
Gaussian kernel in order to suppress noise and match the 
slightly worse g-band seeing, then fitting this smoothed 
i-band LRG image to the g-band image and subtract- 
ing it. The difference images are then smoothed with a 
Gaussian kernel of FWHM « 0'/4 in order to bring out 
coherent features; they show strong evidence of multi- 
ply imaged, relatively blue sources. SDSSJ0037 seems 
likely to be a double-image lens. It is tempting to iden- 
tify quadruple-image morphology in the SDSSJ0216 dif- 
ference image, but the level of noise recommends cau- 
tion. Narrow-band imaging or integral-field spectroscopy 
should permit quantitative study of these systems and of 
those for which broadband imaging is inconclusive. HST 
imaging of any confirmed lenses using the narrow-band 
ramp filter set of the Advanced Camera for Surveys could 
also be pursued to obtain a highly resolved picture of 
lensed line-emitting regions and permit even more de- 
tailed study of the lensing mass distributions. 

The degree of lensing that will be present in a partic- 
ular system will depend on the impact parameter (angu- 
lar offset) of the background galaxy in the source plane 



relative to the center of the LRG: the smaller the im- 
pact parameter, the more lensing will be seen. To ob- 
tain a rough estimate of the number of lenses within 
our sample, we can guess at the unlensed surface bright- 
ness distribution of our background galaxies, compute 
lensed images under an assumed lens model, smear to 
account for seeing, and integrate over the 3"-diameter 
SDSS spectroscopic fiber, then compare to the [O 11] line 
flux values that we have observed. To interpret the re- 
sults in terms of lensing probabilities we must also invoke 
an [O II] AA3728 luminosity function (LF). Appendix 151 
describes the details of a lensing probability calculation 
of this nature that we have carried out. The results sug- 
gest that a total of 19 out of 49 systems are likely to be 
strong lenses — that is, we expect approximately 19 sys- 
tems to have source-galaxy impact parameters less than 
the critical value for multiple imaging. It is important 
to recognize, though, that the "lenses-or-not" question 
does not have as straightforward an answer for extended 
sources as it does for point sources. Different regions 
of our background galaxies will be lensed by different 
amounts, and in general some but not all of the galaxy 
can be multiply imaged. It is more appropriate to ask 
how much lensing is present in any given system. 

There is a sense in which the unknown degree of lens- 
ing in our syst ems may prove to be of scientific interest. 
iBinnevI {2003) raises the concern that current strong- 
lcnsing evidence for massive dark-matter halos in ellip- 
tical galaxies may reflect a selection bias. We quote di- 
rectly: "even if only a minority of ellipticals have massive 
dark halos, nearly all the observed lenses will belong to 
that minority." Our LRG sample, in contrast, has not 
been selected for lensing, but rather for back-lighting. 
The true incidence of lensing within our sample will thus 
constitute a test of whether or not the dark halos in 
currently known lenses are in fact generic to early-type 
galaxies. 

We present a comparison between our lens candidate 
LRGs and the full 51,000 LRG spectroscopic sample 
in Figure El The lens candidate LRGs seem some- 
what skewed towards brighter magnitudes. The con- 
servative interpretation is that the broadband selection 
properties of fainter LRGs are more easily perturbed 
by background galaxies, but it may also reflect a lens- 
ing signal, with more massive galaxies providing more 
magnification. If the logarithmic slope of the under- 
lying [O n] LF is steeper than —1 at the luminosities 
probed by our survey, then magnified lines of sight should 
show a statistically enhanced n u mber of [O i ll emi tters 
ijTurner. Ostriker. OottllT98l . IHogg et a,U lll99l find 
the logarithmic [O n] LF slope to be steeper than —1 
for line luminosities > 10 42 ergs s _1 , whereas the me- 
dian observed [O n] luminosity in our sample is of order 
10 41 ergs s" 1 (again, we have taken Ho — 70 km s _1 
Mpc -1 ). These numbers suggest that magnification bias 
is not the explanation for the observed brightness of our 
lens candidates relative to the full LRG sample, but a 
more definite statement must await follow-up observa- 
tions that capture the total background line flux. 

Most known gravitational lens systems have been se- 
lected on the basis of some combination of source prop- 
erties and lens cross section, whereas the SDSS LRG 
sample is selected based on colors and magnitude. It is 
therefore of interest to compare known early-type lens 
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Fig. 4. — 40" X 40" SDSS r postage-stamp images of LRG systems with confirmed background emission. Gray-scaling is linear from — 3tr 
sky noise (white) to LRG peak brightness (black). Arrows point North; East is 90° counterclockwise from North (i.e., as seen on the sky). 



galaxies to our candidate lenses, although this is dif- 
ficult since we do not currently have truly compara- 
ble observations of the two samples. Figure [7] presents 
our best attempt at such a comparison for lens veloc- 
ity dispersions and apparent magnitudes as a function 
of redshift 10 , with HST known- lens data taken from 
iRusin et al.1 l)2003al) . The known- lens velocity disper- 
sions are estimated fr om lensed image separa tions in the 
manner described bv iKochanek et al] l)2000(k they may 
be systematic overestimates if these lenses are superim- 

10 We plot these quantities as a function of redshift to avoid the 
issue of evolutionary and k corrections. 



posed on the "mass sheet" of a high-density environ- 
ment (IHolder fc Schechterl |2003) . Magnitude compari- 
son is made by transforming SDSS g- and r-band mag- 
nitudes to Johnso n-Morgan F - band using the observed 
transformation of Smi th et alJ l)2002|) . We see that our 
LRG lens candidates are in general of greater velocity 
dispersion than known lenses, and in the redshift range 
where the two samples overlap, the LRGs are more lu- 
minous. The brightness of the LRGs combined with the 
relative faintness of the background galaxies in our sam- 
ple suggests that any confirmed LRG lenses would be 
well suited to the type of detailed lens stell a r-dyna mical 
studies described by iTreu fc Koonmansl (|2002|) and 
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Fig. 4. — Continued 



IKoopmans fc Treul (|2003|) . 

5. CONCLUSIONS 

We have presented a catalog of candidate strong 
galaxy-galaxy gravitational lens systems detected spec- 
troscopically within the Sloan Digital Sky Survey. These 
systems have known foreground and background red- 
shifts; only the detailed spatial alignment of foreground 
and background galaxies remains unknown. We plan to 
conduct integral-field spectrographic follow-up observa- 
tions of these systems, which will allow us to resolve the 
spatial distribution of the background nebular line emis- 
sion that we have detected. These observations should 



confirm a substantial number of lenses within our candi- 
date sample, with many lensed galaxies at lower redshift 
than any other currently known lensed extra-galactic op- 
tical sources. Any lenses confirmed within our sample 
will be of considerable interest for the study of early- 
type galaxy mass distributions, and could have implica- 
tions for lens-time-delay Hq measurements and anoma- 
lous quasar-lens flux ratios. Our sample demonstrates 
the feasibility of the emission-line-based spectroscopic 
lens search technique within the SDSS and other red- 
shift surveys, and we plan to extend the search to higher 
source redshifts in the near future. We have also de- 
veloped and applied a method for abating the influence 
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Fig. 5.— MagIC 8^'8 X 8'.'8 images and difference images of lens candidates SDSSJ0037 (top) and SDSSJ0216 (bottom). Left-hand 
panels show 120-s g-band images of the two systems; right-hand panels show smoothed g — i difference images created as described in the 
text. Note the evidence for multiply imaged, relatively blue objects in the difference images. (In these images, North is left and East is 
up — reversed parity from what is seen on the sky.) 



of night sky emission-line residuals in the 7000-9000-A 
range that allows us to detect [O III] and H/3 emission 
over the redshift range z ~ 0.4-0.8 without an excess of 
false-positive detections. 
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Fig. 6. — Comparison of lens-candidate LRGs (crosses) to full 
LRG sample (gray dots) in the magnitude-redshift plane. The three 
black lines show the redshift-dependent median de Vaucouleurs r 
magnitude of the full sample and the 84th- and 16th-pcrcentile r 
values. 
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Fig. 7. — Comparison of candidate LRG lens galaxies (crosses) 
with known early-type lenses with known lens redshifts (filled 
squares). Top panel: comparison of velocity dispersions er„ mea- 
sured for LRG lens candidate galaxies and inferred from image 
separations A9 in known early-type lens systems wit h known lens 
and sou rce redshifts. Image separations are taken from lRusin et ahl 
1 2003a), and conv ersion to velocity dispe rsion estimates is made 
via Equation 5 of Kochan ck et alJ 120001) . For visual clar i ty, we 
omit error bars on the latter data points; [Kochanck et al.j J200Q]) 
estimate 10% errors in the A9 — * o v conversion. Bottom panel. V- 
band magnitude comparison. SDSS g- and r-band magnitudes are 
c onverted to Johnson -Morgan V via the observed transformation 
of lSmith et~al"l f2002'l: we have also correct ed for an 0.2-magnitud e 
systematic error in the SDSS magnitudes (Abazaiian et al .1 12003(1 . 
Error bars on SDSS photometric quantities are smaller than 
plotted sym b ol size s. Known-lens photometric data taken from 
IRusin et all I2003al) . All magnitude s have been correcte d for 
galactic extinction using Schlcecl. Finkbcincr. & Davis (1998) dust 
maps. 



22 



Bolton et al. 



APPENDIX 
A. NOISE MODELING 

If the model of a purely Gaussian noise spectrum described by ax were correct, then the distribution of scaled 
residual specific fluxes 

xx = f^/ax (Al) 

across all spectra would be Gaussian with unit variance for all wavelengths A. This is unfortunately not the case in our 
sample. Imperfect night-sky emission-line subtraction and other miscellaneous effects give rise to an excess of "high- 
significance" outliers beyond the predictions of a Gaussian model, leading to a deluge of false-positive astronomical 
emission-line candidates when the procedure described in § 13. H is applied, particularly in the 7000-9000-A region of the 
spectrum where the [O in] A5008 line at redshifts z w 0.4-0.8 appears. The most drastic solution is simply to mask all 
sky-afflicted wavelengths. Rather than concede such vast spectral coverage (which would drastically reduce our survey 
volume) , we describe the observed distribution of scaled residual specific fluxes xx within the LRG sample with a more 
detailed empirical noise model. The generally Gaussian behavior of scaled residuals at low significance combined with 
the excess of high-significance residuals is well described by a mixture of Gaussian and Laplace distributions, expressed 
parametrically as 

p(x) dx — [aexp(-a; 2 /2f7g) + bexp(-|ar|/o- e )] dx (A2) 

(For history and applications of the Laplace distribution, see lKotz. Kozubowski. fc Pode orski 200l]) The parameters 
of this distribution are wavelength-dependent, but we suppress this dependence in out notation. The values of a and 
b are related by normalization: 

/+oo 
p(x) dx = V2na g a + 2<r e b = 1 . (A3) 
-oo 

We also fix the following relations between parameters, based on strong correlations observed in free-parameter fits to 
the distribution at each wavelength: 

a e = a g -0.38 , (A4) 
6 = 0.09 x o g (a + 6) . (A5) 

The result is a one-parameter noise model to fit to the distribution of xx across the sample at each wavelength. (The 
numerical values 0.38 and 0.09 are fixed by minimizing the sum of binned \ 2 values for fits across all wavelengths.) We 
relax conditions (|A4|) and (|A5|) and fit freely for a e and b at a few isolated location s in the s pectr um, where the effects 
of sky-subtraction residuals are especially strong and the correlations that suggest (|A4|) and (|A5|) break down — regions 
near 5577 A, 5894 A, 6305 A, and 6366 A. Additionally, some regions of some spectra are characterized by extreme 
and correlated excess variance, so for each spectrum we convolve \xx\ capped at 5 (to limit the influence of single 
pixels) with a 100-pixel boxcar filter and exclude from the noise-modeling sample any pixels within a boxcar whose 
value exceeds 1.25. 

We use our fitted noise model to re-scale the reported ax values such that the new distribution p{x) dx of scaled 
residual flux values at each wavelength becomes Gaussian, while preserving the position of individual x-values within 
the cumulative distribution, then proceed as described in S I3.1I Both the reported noise ax and the measured residual 

(r) 

flux values contain information about the actual error in the presence of imperfect subtraction, so it is sensible 
to base an effective noise rescaling on their ratio xx in this manner. By fitting the noise distribution parameters 
independently at each wavelength, we also capture the localized effects of individual night-sky lines. 

B. LENSING PROBABILITIES 

Our strategy for assessing lensing probabilities in our sample centers on the construction of an approximated prob- 
ability density p(b) db for the unknown impact parameter b of the background galaxy in each system. The following 
observed quantities are input to the calculation: the LRG and background redshifts, the LRG velocity dispersion, 
the background [O n] AA3728 line flux received by the 3"-diameter spectroscopic fiber, and the median seeing for the 
spectroscop i c plat e under consideration. We also make use of the [O n] line luminosity function (LF) reported by 
iHogg et all (|1998T> . We adopt the same SIS model for the LRG mass distribution as was used to obtain the AO values 
in Tabled an d we model the background galaxies as exponential disks with a half-light radius of approximately 3 kpc 
(fixed to 0'.'5 at z = 0.5). 

For each system, we explore a range of impact parameters b from to 5". At each 6-value, we generate a lensed 
image of the model background galaxy, then convolve it with a Gaussian point-spread function corresponding to the 
median seeing reported in the spectroscopic plate header. We then integrate the image over a 3"-diameter circular fiber 
aperture centered on the model lens. The result is a tabulated function f(b) giving the fraction of the intrinsic flux 
received by the fiber; that is, if the total galaxy line flux were S in the absence of lensing and limited fiber sampling, 
the [O n] line flux received by the spectroscopic fiber from a background galaxy with a source-plane offset b would be 
Sab = f {b)S. In general f(b) may be greater or less than one due to the competing effects of lens magnification and 
incomplete sampling by the fiber. 



SDSS Spectro Lens Search. I. 



23 



Next we adopt the [O n]-emitter LF reported bv lHogg et alJ l)1998(l by fitting a Schechter function to their Figure 6. 
After converting from logarithmic units, the number of [O II] emitters per unit volume in an interval dL at line 
luminosity L is well approximated by 

4>{L) dL (x L a exp(— L/L*) dL , (Bl) 

with a ~ —1.3 and ~ 3.4 x 10 ergs s _1 (the overall normalization is unimportant for our purposes). We 
make a crude conversion from their assumed (Ojifjfl/y) = (0.3,0) universe to our cosmology by scaling their reported 
luminosities up by a factor of 1.2: the ratio of squared luminosity distances in our cosmology to theirs ranges from 
1.16 at z = 0.3 to 1.27 at z = 1, and the bulk of galaxies in their study fall within this range. Assuming the form of 
the LF does not evolve, it corresponds to an intrinsic flux function at any redshift z for the number of galaxies per 
unit redshift per unit solid angle within some intrinsic flux range dS about S: 

^(z,S)dS = N(S,z)S a exp(-S/S±)dS , (B2) 

with the same a as the LF and S* = L^/[47r£)| / (z)], where Dl(z) is the luminosity distance to redshift z. N(S, z) is a 
flux- and redshift-dependent normalization, the form of which will prove unimportant. 

We can now derive a joint probability density function (PDF) for the observation of an [O ii]-emitting galaxy behind 
a given LRG at impact parameter b, redshift z, and with line flux Sfib in the fiber by setting the differential probability 
proportional to the corresponding expected differential number count and making use of the known relationship of S^b 
to intrinsic flux S through /(&): 

p(b, z, Sah) db dz dSfib ocip(z,S) dft dz dS (B3) 
= i>[z, 5 fib //(6)] (2tt6 db) dz [dS &h /f(b)} . 

The term dfl = 2nb db represents the solid angle in the source plane of an annulus of radius b and thickness db. The 
observed quantities S D bs and z b s for the system are equal to the system's true Sfib- and z-values plus some observational 
noise that is independent of b, so assuming the noise is small relative to the scale on which the joint PDF varies, we 
may reinterpret the joint PDF as an approximate conditional PDF on b given z bs and S'obs: 

p(b; z ohs , S'obs) db = N'(z ohs , S obs ) b [/(b)]" 1 ip[z ohs , S ohs /f(b)} db . (B4) 

The normalization ^'(zobs, S b s ) need not be derived explicitly, since we can simply compute the right-hand side 
without it for the relevant range of 6-values and normalize afterwards. With this PDF in hand, we can finally assign a 
"lensing probability" to the system as the integrated probability for impact parameters less than the critical value for 
multiple imaging in the SIS model — the "Einstein radius", equal to one-half the A6* value given in Tabled Summing 
this probability over all systems gives the estimate quoted in § 14.21 of 19 strong lenses out of 49 candidates. 
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